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Abstract 

Using a statistical approach in the framework of non-covariant perturbation 
theory the distributions for light and charmed quarks in the hadron have been 
derived, taking into account the mass of the charmed quark. The parameters of the 
model have been extracted from the comparison with NA3 data on hadroproduction 
of J/tp particles. A reanalysis of the EMC data on charm production in muon- 
nucleon scattering has been performed. It has been found in comparison with the 
conventional source of charmed quarks from photon- gluon fusion, that the EMC data 
indicate the presence of an additional contribution from deep-inelastic scattering on 
charmed quarks at large x. The resulting admixture of the Fock states, containing 
charmed quarks in the decomposition of the proton wave function is of the order 
of 1%. The approach presented for the excitation of the Fock states with charmed 
quarks can also be applied to states with beauty quarks as well as to the hadronic 
component of the virtual photon (resolved photon component). 



1 Introduction 

The production of heavy flavours in lepton-hadron and hadron-hadron collisions is a very 
important tool for a quantitative test of QCD and for searches for new physics. Due to 
the presence of the point-like probe particle (lepton) and the possibility to control the 
QCD scale of the hard sub-process deep inelastic scattering has a number of advantages in 
comparison with hadronic reactions in the analysis of charm production. The QCD-based 
parton model has been remarkably successful in describing a wide variety of high-energy 
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processes involving energy scales much larger than the masses of the known particles 
and the partons themselves. Many analysis of charm production performed within the 
framework of the parton model assume that the hadron consists only of massless or 
approximately massless partons (gluons, u, d and s quarks). The heavy quarks (charm, 
bottom) are treated as massive objects which are external to the hadron. In the DIS- 
neutral current reactions this kind of consideration leads naturally to the 0(a]) "photon- 
gluon fusion" (PGF) mechanism 7$ — > cc, as the main mechanism for heavy quark 
production. In hadronic collisions the analogue "parton-parton fusion" processes, gg — > 
cc, qq — > cc, are expected to contribute. These parton fusion processes are flavour 
creating (FC), since the heavy flavour is created by the interaction with a light constituent 
of the hadron. 

The existing experimental data on fip collisions jl| show some irregularities which 
are inconsistent with the PGF predictions. The experimental observation of a deviation, 
of the charm distribution at large pseudorapidities (i.e. charm production close to the 
direction of the proton beam) from the conventional predictions in ep scattering have 
also been reported by ZEUS experiment at HERA collider ||. In hadronic collisions 
the interpretation of the data on open charm production within the standard parton- 
parton fusion scheme followed by hadronization of charmed quarks meets problems for 
charmed particles at large xp. In this region the charm distributions are harder than then 
predictions of the factorization approach. Furthermore, the yield of charmed particles 
containing the valence quarks of the projectile significantly exceeds the yield of their 
anti-particles. Models considering recombination of the newly created charmed quark 
with one of the valence quarks of the projectile [§, || or string fragmentation || can 
improve the situation with open charm production. In both approaches a part of the 
proton remnant momentum is imparted to the final state charmed particle increasing its 
momentum and improving the agreement with the experimental observations. However 
these models have problems with describing J/ip and double- J/ip production as well as 
the A-dependence of the charm production cross section in hadron-nucleus collisions at 
large x F |J, | . 

A part of the discrepancies between data and models can be resolved by introducing 
the scheme of flavour excitation (FE), which assumes that also heavy quarks can be 
constituents of the hadron. We note that considering heavy quarks as external to the 
hadron is appropriate when the characteristic scale of the process (/1) is less or of the 
order of the mass of the heavy quark, i.e., \x < tuq. This condition holds for c and b 
quarks for most fixed-target experiments. The ep collider HERA gives an opportunity 
to investigate the heavy quark production at Q 2 scales much larger than 4mg. At such 
scales it seems justified to consider charm (and bottom) quarks as light objects. 



The so called "variable flavour number scheme", suggested in papers R ITDL O 



com- 



bines both, the FC and FE schemes and ensures a "soft" transition between the two 
production mechanisms. A more sophisticated approach for the FC mechanism re-sums 
effectively the large logarithms of the type [a s (/x) ln(/i 2 /mg)] n which limit the validity 



2 



of the conventional FC mechanism to the region \x ~ O(ttiq). It has been shown in 
these papers that the contribution from the scattering on a constituent charm quark of 
the proton (process of order 0(az®)) becomes more important than PGF at x > 0.1 al- 
ready at surprisingly low Q 2 > 20 — 30 GeV 2 . The authors used splitting functions and 
standard distributions for massless partons, including heavy flavours and introduced the 
" slow-rescaling" variable x — > x [1 + (mq/Q) 2 ] to account for the charmed quark mass. 
However perturbative QCD requires the scale fi to be large and needs as input the initial 
parton distributions for their evolution. The initial distributions for heavy quarks are 
not necessarily similar to the distributions of light partons due to non-zero quark mass, 
which is comparable with the QCD scale. This kind of consideration is closely related to 
the old question: which type of high order QCD corrections have to be assigned to the 
matrix element and which due to QCD evolution of the parton distributions. So far we 
have no clear understanding of this problem. 



The authors of |L2| have suggested a procedure to obtain the distribution of massive 
charmed quarks in the hadron. They considered the proton wave function decomposition 
which may contain the | uudcc > Fock state component called "intrinsic charm" (IC). 
Such a state may appear as a quantum fluctuation of the hadron wave function and may 
become free in interactions with substantial momentum transfer. In this case the proton is 
described as a decomposition in terms of colour-singlet eigenstates of the free Hamiltonian 
\uud >, \uudg >, \uudqq >, .... Over a sufficiently short time the proton can contain 
Fock states of arbitrary complexity, including pairs of charmed quarks. In the proton rest 
frame the life-time r of such fluctuations is of the order of the nuclear time ~ Rh, where 
Rh is the hadronic size. On average there are in addition to the valence quarks extra 
partons (gluons and qq pairs). In the infinite momentum frame a partonic fluctuation 
will be "frozen" and can be observed, e.g., in lepton-hadron scattering. The charmed 
quarks are heavy objects and their life-time is much smaller than for light partons. Thus, 
on average, the admixture of heavy quark pairs is expected to be small, ~ (m q /mq) 2 . 
Because the quantum fluctuations in the initial proton is defined by the colour field self- 
interaction, the structure of Fock states of the proton can be considered independently of 
the hard interaction, providing the initial non-perturbative parton distributions. These 
distributions will evolve in hadron- or lepton-hadron collisions due to large momentum 
transfer. 

In the present paper we modify and generalize the statistical approach to the Fock state 



hadron structure with heavy quarks, suggested in Ref . |12], [13], [14] , within the framework of 
the non-covariant perturbation theory. We obtain scaling expressions for the heavy and 
light parton distributions in the infinite momentum frame. We calculate the charmed 
structure function of the proton, F2°\x,Q 2 ), taking into account the QCD radiative 
corrections order of ai as well as the mass corrections caused by the non-zero values 
of the c quark and the proton masses. We use the experimental data for it A — ► J/tpX 
and p,p — > /iccX H to evaluate the parameters of the model and present the relative 
contributions of the PGF and the IC mechanisms to the charm structure function of the 



3 



proton. Note that we use the terms FE and intrinsic charm (IC) for the charm production 
mechanisms involving constituent charmed quarks of the proton. 



2 Description of the model 

2.1 General features 

In QCD high energy hadrons are coherent superpositions (Fock state vectors) of quarks 
and gluons. Note that the lifetime of the fluctuation At ~ 1/AE m 2P h / (M 2 — m 2 ) 
(Ph is the momentum of hadron, m is its mass and M is the mass of the fluctuation) can 
be very large at high energies even for large mass values of the fluctuation. 

Based on above picture of the proton and using the statistical approach, the authors of 
Ref . [O] achieved a good description of the proton structure function. In Ref . flTl] we have 



presented a statistical consideration of the hadron structure and obtained non-invariant, 
i.e. frame dependent, expressions for parton distributions. In principle the frame depen- 
dence can take place at sufficiently low energies while in the infinite momentum frame 
one expects invariant expressions. 

We take all partons on the mass shell and use non-covariant perturbation theory. 
Thus we consider a hadron as a statistical system which consists of N quarks carrying 
quantum numbers of a hadron, two charmed quarks c, c and a system of n light partons 
(gluons and quarks) carrying in total the quantum number of the vacuum. 

In non-covariant perturbation theory the probability to produce a m particle final 



state in the case of the instant interaction potential looks as follows 15]: 



where c?$^ is the element of the m particle phase space, m = N + 2 + n, Ph and Eh are 
the momentum and energy of the considered hadron; Pf in and Ef in are the momentum 
and energy of the final state partonic fluctuation, respectively; d&J^ describes the Lorentz 
invariant phase space: 



= n^. (2) 

i=i 

The 5-function ensures conservation of the total 3-momentum. 

Due to the sharp cut-off for the transverse momenta of partons it is sufficient to 
consider only the longitudinal phase space: 
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d 3 p <g . . 

where £ = p z /Ph, p = m±/Ph and m± is the transverse mass of the parton. 

Following the parton model we assume independent primordial distribution of each 
parton and replace: 



m 
i=l 

where is the probability density to observe the parton with the momentum fraction 
£. Therefore, the probability to observe an m parton Fock state looks as: 



p\ m / m \ 



(5) 



The omitted common factors will be incorporated in the general normalization. The 
integration over m parton momenta can be performed |13| with help of the integral rep- 
resentation of the ^-function: 



-oo 



2tt5(x) = / due wx . (6) 

J — oo 

If we perform the integration over all £j, we obtain: 

i r+oo m 

= — / due iv Y[piM. (7) 

Z7T »=i 

Here we extended integration of £ to oo due to the presence of the fast oscillating exponent 
and defined as p%{v) the Fourier transformation of the parton density Pi(£) : 



••DC' 



'0 



Following [HI we introduce different probability densities for valence quarks, charmed 



quarks, gluons and light sea quarks, p v , p c , p g and p q , respectively Because all sea light 
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partons of the same type have the same distributions we have to sum over all possible 
permutations of the n light sea partons (gluons and qq pairs separately). Thus, we have 
for the probability to observe a Fock state with N valence quarks, one cc pair, n g gluons 
and n q pairs of light sea quarks (n = n g + 2n q ) the expression: 



w 



(n) 



N 



= 2tt / dv e w [p v (v)] N pcciy 



E 

n g +2n q =n 



p 9 (yY B Pqiy) 



2n n 



n, 



(2n q 



(8) 



where p c5 is the probability to create cc pair. 

The factors l/n g \ and l/(2n q )\ take into account the nondustinguishability of gluons 
and quarks, respectively. To obtain the total probability for all Fock states, containing 
the cc pair one needs to perform a summation over < n < oo. This summation can be 
carried out by using the properties of binomial sums: 



E 



rig 2n q 
Pg Pq 



' n a \ (2n a )\ 

n g +2n q =n 9 V q} 



1 

2 



P- 



ni 



+ 



pi 



p± = p g ± p q 



(9) 



The general form of the statistical sum is [16]: 



4 C) = £Mr } ( 10 ) 

n=0 




dve w \p v {v)] N p C c{v) exp[p 9 (z/)] cosh[p q (u)]. 



It is clear from Eq.(|5|), that the distribution Pi(£) for the ith parton can be obtained 
if one omits the integration over momentum of the considered parton. In general, one- or 
many-particle distributions can be derived from the statistical sum taking the functional 
derivative of the desired function(s) fl6f . Thus, the inclusive distribution of light partons 
and the distribution of the cc pair, normalized to unity, looks as follows: 



Sp, 

Pacific: Cc) = ~Z(cj Pcci^c: Cc) JTT = T^cf Pcci^ci Cc) Ccc(l — Cc — Cc; - 

£ at " cc & i\r 



r{o) 



(11) 



where, Cj(l — £) and C C c(l — £ c — £5) are the correlation functions, ensuring the momentum 
conservation. 
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2.2 Probability densities and parton distributions. 



The origin of the cc pair is the same as of light sea quark pairs, namely the splitting of 
the gluon into a virtual cc pair, g — > cc. In principle, at sufficiently small energies, when 
icPh < m ±c, one can expect p c ~ const, but in the limit Ph — > oo, one has case £ ^> u c 
and therefore a situation similar to light sea quarks. We shall compare the model with 
fixed target experimental data, neglecting the transverse mass of the charmed quarks for 
the projectile hadron. Thus, following |HJ we can represent the probability densities 
as: 



Pc(0 OC 

for valence quarks, gluons, light sea quarks and charmed quarks, respectively, with a g and 
a q being unknown constants. 

In the infinite momentum frame we neglected the transverse mass in the probability 
densities of the valence and charmed quarks, Eq . (|i~2l) . At the same time for sea par- 
tons we hold temporarily the term \x in the denominators to perform later the Fourier 
transformations. In the final expressions we shall put fi — > 0. 

From experiment we know that the valence quark momentum distribution at small 
£ is approximately proportional to l/v 7 ?, i-e., a = 0.5. We shall use this value in our 
comparison with experimental data, but in the formulae we use the general expression 

AD- 

Let us consider the energy denominator in Eq.(|I|). Taking into account the momentum 
conservation Ph = ^\ pi, from the light-cone expansion in the infinite momentum frame 
one obtains: 



2 ' 
9 



1 



2 ' 



i 



(12) 



E 



fin 



E h 




E 

i=i 



m 



Li 
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(13) 



Because the charmed quark transverse mass is much larger than Mh (the mass of the 
hadron) and than the transverse mass of the light partons m±i, the energy denominator 



will be proportional to [E| : 
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oc 



ScSc 



(Efin — E in ) (£c + £c) 2 



(14) 



The validity of the approximation for the light-cone expansion of the energy denomi- 
nator for Fock states with heavy quarks has been considered in Ref . |14[ . 

We can introduce the expression (14) in the definition of the probability density p c5 
to observe the pair cc, 



2c2 



tec 
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Similar to case of the valence quarks we keep the more general form for cc probability 
density in further formulae: 



We introduced in above formula the phenomenological parameter (3 to take into ac- 
count possible deviation of the charm distribution from the phase space approximation 
Eq.(^) at moderate energies or momentum transfers. In analytical expressions we use 
generalized formula Eq . (|15D . But in numerical calculations and in comparison with data 
in present paper we use the phase space approximation and hold (3=1. 

Therefore from Eq.fllOD we obtain for the statistical sums for Fock states with (Z^) 
and without (Zn) the cc pair in the limit fi — > 0: 



where T(x) is the gamma function and g = a g + a q is the unknown parameter of the 
model, characterizing the level of the sea in the hadron considered. 

Note here, that the integral for light sea partons is ~ and logarithmically diverges 
in the limit \i — > 0. But this divergence can be incorporated in the general normalization, 
as seen from Eq . (|TTj) and is not important. 

The analytical expressions for the statistical sums with and without the cc pair are: 



N 



N 



j_ r(2/3) f (J3) 

v g ^ 5 



1/9 ' 



(16) 
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(c) = [T(a)] N rm m 
N T(aN + 2/3 + g) 

Z - [T{a)]N (17) 

r /2 (sin 99 cos 

/cW = / ^ 

Jo (sin + cos 

For integer values of (3 the function f c ((3) is: / c (0) = 1, / c (l) = 1/6, / c (2) = 1/30 and 
/ c (3) = 1/140. For arbitrary values of /3 the integration can be performed numerically. 

The result for the parton momentum distributions for the valence quarks, sea partons 
and charmed quarks is: 



vis = % e- 1 (i - o a{N - 1)+2P+g - 1 - 



z N 



S(0 = 9 T 1 (1 - Q ^ 2 **- 1 , 

P cc -fe&) = (i - £ c - , (18) 

C (0 = ^^(i-o aAr+ ^4 c) (0, 

^ - Jo <*!/ • 

The valence and c quark distributions are normalized to unity P(£) = 1. 

Despite the fact that the mass of charmed quarks does not enter directly the final 
expression for the probability of the Fock state, we see from Eq . (|i"8|) that due to the factor 
£p with (3 > the distribution of charmed quark is much harder than a distribution of 
light sea parton. The origin of this hardness is the large value of the charmed quark mass. 
This effect must be taken into account by any phenomenological parameterization of the 
initial (not QCD evolved) charmed quark distribution. 

If the SU(3) symmetry in the sea is broken one can introduce a suppression factor X s 
for strange quarks and obtain: 
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(19) 

Since the total probability for the excitation of a cc pair in the hadron is unknown we 
assumed the normalization to be unity. The final result can be obtained by multiplication 
with the factor N c , which can be extracted from the experimental data. In the next 
sections the model is compared with the experimental data in order to evaluate the free 
parameters N c and g. 

Let us make a comment about the excitation of the states with beauty quarks. If we 
consider b quarks within the framework of this approach, we obtain again the expressions 



(15) and (0) neglecting the mass of the charmed quark in comparison with the much 
larger mass of the beauty quark. In this case, the charmed quarks have to be treated as 
massless partons of the sea. 

One can apply this model also to the excitation of the heavy flavours in the virtual 
photon ("resolved photon") in e + e~ annihilation or in photoproduction. To obtain the 
necessary distributions it is enough to omit the term \p v {v)\ N in Eq . (|T0l) because, obvi- 
ously, the photon does not contain valence quarks. 



3 Comparison of the model with J/ip hadroproduc- 
tion 

It was noted in the introduction, that there exist a number of hadronization mechanisms 
which more or less successfully describe the inclusive open charm yield at large xf- These 
schemes incorporate some aspects of the hadronization process, but there is no commonly 
accepted mechanism based on a well founded theoretical approach, which describes all 
available data. The uncertainties in the existing hadronization models are too large to 
perform a direct evaluation of the parameters of the model for charm quark distribution 
presented here from the experimental data on open charm production. 

On the other hand the experimental data from NA3 J7| on the production of J/ip par- 
ticles in hadron-platinum collisions indicate an unusual production mechanism. The au- 
thors of NA3 identified two different components in the xp distribution of J/tp mesons: a 
"hard" component with the usual A-dependence a a ~ Acr^ and a " diffractive" com- 
ponent with a weaker A-dependence, namely, <ja ~ A 0,77 ax for incident pions and 
a a ~ A°' 71 <jn for incident protons. The relative contribution of the "diffractive" com- 
ponent is ~ 0.20 and ~ 0.30 for pion and proton projectiles, respectively. The "hard" 
component, as shown in the same paper 0, can be described well by the conventional 
QCD parton-parton fusion mechanism. The linear A-dependence agrees well with pre- 
dictions of the QCD factorization theorem valid for hard processes. 
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The unusual A-dependence of the " diffractive" component can naturally be described 
in the model of intrinsic charm [[R] and the model pi based on the Gribov approach 
to particle interactions with nuclei. It can be shown || that for J/ijj production at NA3 
energies, P[ a b ~ (150 — 300) GeV, a deviation from the A 1 behaviour for inclusive spectra 
only for those components of the initial partonic configuration, which contain the heavy 
state. Therefore, we can ascribe the "diffractive" component in J/ip distribution, seen in 
NA3 experiment, to the interaction of the hadronic Fock state, containing a cc pair and 
to use these data to evaluate the parameters of our model. 

To estimate the longitudinal distribution of J/ip particles we use the recombination 
In this model the differential cross section for J / ip production can be written 



model 

as: 
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= CFtot [ d^d^Sce^ctc) R(£c,&iXf) , (20) 



dx 



where, S c5 (!; c , £g) is the two-particle distribution of c and c quarks with momentum frac- 
tions x c and x s , respectively; i?(£ c , £ 5 ; Xi?) is the recombination function describing the 
probability for two quarks to coalesce in the final J/ip meson with momentum fraction 
xf- In the simplest case 

R{£c,&xf) = 5{x F - £ c - 

ensuring longitudinal momentum conservation. In principle, for charmed particles, the 
primordial transverse momenta of the initial c quarks can reach large values (> 1 GeV) 
and have to be taken into account (see, e.g., |§]). For the aim of the present paper 
this is not important. Since we consider only longitudinal distributions the integration 
over transverse momenta is included in the total normalization on experimental data. 
As a result, from Eqs.fllBD, (p0|) one obtains for the xp distribution of J/ip particles the 
following expression: 



dcr(J/i>) ™ Z N 1)] m_! , aN +g-l /,,-, 

~H7~ M 4) W& + 2) (1 " Xf) ■ (21) 

For the fit we used only 7r iV data from assuming for the number of valence quarks 
N = 2. We do not use the pp — > J/ip data because we actually don't know the probability 
for the cc pair to form a J/ip meson. This probability can be different for incoming 



7r-mesons and protons. Therefore, expression (21) has two free parameters: the total 
normalization and the parameter g characterizing the level of sea partons in the 7r meson. 
The fit gave the following result: 
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g = 1.35 ± 0.09, 



(22) 



with x 2 /NDF « 0.9. It should be kept in mind, that the parameters a and (5 have been 
fixed, a = 0.5, /3 = 1. Results are shown in Fig.|I](a-d). 

To recalculate the predictions of our model to NA3 data on pA — > J/ip we have 
to put N = 3 in Eq.([H]) and to change the normalization o^f, found for the ir data, 
by the W pn a^[, where W p7T = atatipp J W) / 'vtoti^p — * J/ip)- To estimate the ratio 
W p7T , we neglect any dynamical effects and use the combinatorial consideration based on 
the following simple assumptions. We consider fast quarks only, i.e., the valence and c 
quarks of the beam particle. Each massless quark has two spin states and three colour 
states (we neglect masses for fast c quarks). To form white, spin 1 final J/ip meson, we 
have to choose c ans c quarks with parallel spins in colourless state. It is clear that the 
probability of random choice of the cc pair with necessary quantum numbers is being 
given by the binomial coefficients G% , where n is the total number of states in the beam 
particle. The statistical weight of the final state is the same for both beam particles and 
cancels in the ratio. In 7r-mesons we have four fast quarks (udcc) and, consequently, 
n n = 3x2x4 = 24 states. In the proton we have five fast quarks (uudcc) and 
n p = 3x2x5 = 30 states. Thus, W pn = Cf/Cf = 92/145 in good agreement with 
data. 

In Fig .|^(e) we plotted the resulting da /dxp(J ftp) for pp interactions. As one can see 
from Fig.[TJ(e) , the model gives a satisfactory description also for J/ip production in pp 
collisions. This allows us to use the same parameter g for the analysis of EMC data on 
charm production in muon-proton scattering. 

Let us remark on some irregularities, seen in J/ip distributions in Figs.|l|(a-c) around 
the xp ~ 0.8 — 0.9. E.g., in Fig.|T|(c) the point at Xf ~ 0.85 is about one order of 
magnitude higher, than the theoretical curve. If it is not a statistical fluctuation in data, 
this discrepancy can be easily understood within the framework of the model 0, which 
predicts that at very large A and xf — > 1 the production cross section can have an A 1 / 3 
dependence, if there is a final state interaction of the J/tp particle with nuclear matter. 
In this case we would obtain an additional factor 195 1 / 3 ~ 5.8 in good agreement with 
Fig.|(c). 

In Fig.|I](f ) we also present the distributions x q(x) for charmed and valence quarks in 
the model and the distribution of the valence u quark from the MRS (G) parameteriza- 
tion. 

4 Charm Electro-Production 

4.1 IC Structure Function and Sub— Leading Corrections 

The cross section for charm production in deep inelastic muon-proton scattering is: 



12 



da 2tt a 2 [1 + (1 - yf] 



F 2 c (x, Q 2 ) , (23) 



dxdQ 2 xQ l 



where, x = 2 (Pq) ^ s ^ ne Bjorken variable and y = ^ is the fraction of the 
muon momentum carried by the virtual photon (we neglected the contri- 
bution from the longitudinal structure function F£)- 

The charm structure function of the proton in the present approach can 
be represented as a sum of two terms: 

F<(x,Q 2 ) = F!f GF \x,Q 2 ) + N c Fi IC \x,Q 2 ). (24) 

The first term describes the conventional photon-gluon fusion mechanism [ig —>■ ficc 
(FigJJa) and the second term represents the direct scattering of the muon on the con- 
stituent charmed quark of the proton, fi + c — > fi + c (Fig.|2|b). N c is an unknown 
normalization constant to be found from the comparison with the experimental data. 

Within the framework of the naive parton model the IC structure function is connected 
with the distribution c(x, Q 2 ) of the charmed quark in the proton as: 

F 2 c (x,Q 2 ) = 2elxc( Xj Q 2 ), (25) 

c(x, Q 2 ) is the momentum distribution of the c quark and e c = | is the electric charge of 
the c quark. 



The PGF charm structure function is 18 



Jy/T+AXx K K 



(26) 



where, 



f2(z,Q 2 ) = !^£K Z lv e [-\ + 2z(l-z)(2-\)] 

+ [l-2z(l-z) + 4Xz(l-3z) - 8AV]lnf±^j. 



(27) 
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In the above expressions s = Q 2 (1 — z)/z, V c (s) = y/l — 4m^/s is the velocity of c 
quark in the (7*7) CM system. 

To compare the model with the experimental data we need to take into account the 
dependence of the c(x) on the momentum transfer Q 2 . This dependence comes from two 
sources [IS]. The first one results from the non-zero masses of the proton and the c quark. 
The second source are the first-order QCD radiative corrections (Fig.^|). 

To take into account effects of the non-zero masses we replaced the variable x by the 



variable ( [O 



x - CM = , Vr t^lV (28) 
1 + V 1 + 4px 2 



and the distribution c(x) by the function 



c(x) -> c(C,C) 



The parameters p, x, A, ( are then given by: 




< C < C 



(29) 



P = 



in 



2 j 



A 



1 



l+4A-p' 



)2 ) 



C = C(*> 



(30) 



The first-order correction to the structure function F% (x, Q 2 ) can be represented as 
a convolution of the c quark distribution c((, Q with the radiative corrections. As a result 



the IC structure function with the radiative corrections has the form 19 : 



F 2 c (x, Q 2 ) = 2 el C c(C, C) + 2 e 2 c ( J ^ c(y, ^, (31) 

The expression for the first-order radiative corrections a^\z, A) is given in the Ap- 
pendix. 

In this paper we don't consider the QCD evolution of the constituent charmed quark 
distributions and assume that the QCD radiative corrections to the matrix element of the 
virtual photon absorption give the correct description of the as effects. At large Q 2 this 
point needs a more careful study. There is also the problem of a proper description of the 
QCD evolution of the heavy quark distributions at not too large momentum transfers, 
where the mass of the heavy partons cannot be neglected. 
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4.2 Comparison with EMC data 



In order to evaluate the contribution to the charm part _F 2 C of the proton structure function 
F 2 from the scattering on constituent charm quarks (also called intrinsic charm quarks) we 
used the data on charm production obtained by the EMC collaboration in fip collisions 
at = 200 GeV. The EMC collaboration presented data on a(j*p — > cc + X). Thus, 
we need to extract F£(x, Q 2 ) from the data taking into account the difference in definitions 
of the virtual photons flux, being used in Eq.fl2"B|) and by EMC collaboration [2I| as well 
as the finite size of the experimental bins in (x, Q 2 ) plane. 

According to the Equivalent Photon Approximation (EPA) [F22| the cross section for 
muon-proton scattering can be represented as: 



da(fip —> ccX) = cr(7*_p — > ccX) dn 7 , (32) 

where c?n 7 is the differential flux of equivalent photons. The definition of the equivalent 
photon flux is, to some extent, arbitrary. The conventional form, used in Eq.(^) for 
Q 2 /E 2 < 1, is of the form: 



*i, = £[l + (l-,,)lfg. (33) 



The EMC Collaboration used a slightly different definition for the photon flux [21 
which includes an additional factor (1 — x) in the right-hand-side of Eq.fl33|). Taking into 
account this factor and approximating the differential flux dn^/dxdQ 2 by An/ AxAQ 2 
(where (Az/, AQ 2 ) is the experimental bin), we can relate the structure function to the 
experimentally measured 7p cross section, <J^p P , for charm production: 



F ^ Q ^ = 4^1 + (!-„)»] Q2) (34) 

The value of An 1 is found by integrating the expression for the photon flux over each 
experimental bin (Au, AQ 2 ). Defining Av = v^ — Vx and AQ 2 = Q\ — Q 2 ., we obtain: 



An 7 (z/, Q 



2\ _ 



°L J In ( Q% 



21n^-|( 



V2 ~ V\) + 



2E 2 



(35) 



Q\ Qi 



21n^ + ^ 

v\ hi 



15 



Mass corrections for IC structure functions have been used in the form of ^-scaling 
The radiative corrections (see Appendix) have also been taken into account. The 
strong coupling constant a s (Q 2 ) and the Aqcd has been chosen according to the PDFLIB 
parameterization E3I . 



To perform a fit to EMC data we used expression (33) for the charmed structure 
function with two free parameters: N c being the normalization of F 2 and m c the 
mass of charmed quark, entering the PGF structure function and the radiative and mass 
corrections. For the gluon distribution we have chosen the MRS(G) parameterization, 



which is the default one for PDFLIB 7.09 23 



The results of our fit to the EMC data are shown in Fig.|4j(a,b). The fit yields for the 
admixture from scattering on the constituent charmed quark 



N c = (0.9 ± 0.2) % 

and a value of the mass of 



m c = (1.43 ± 0.01) GeV 

for the charmed quark mass. We have used also other PDFLIB parameterizations and, 
within errors, obtained similar values. 

Figure f|(c) presents the ratio _F 2 /C /F^ GF versus x for some values of Q 2 . We see 
that at large x (x > 0.1) an 1% IC component dominates the charm production for 
Q 2 < (10 - 12) GeV 2 . 



5 Conclusion 



In this paper we modified and generalized the statistical approach to the Fock state 
hadron structure with heavy quarks, using the framework of the non-covariant pertur- 
bation theory. We obtained scaling expressions for heavy and light parton distributions 
in the infinite momentum frame. We calculated the charmed structure function of the 
proton F2°\x, Q 2 ) taking into account the QCD radiative corrections of order as as well 
as the mass corrections to the structure function caused by the non-zero values of the 
c quark and the proton masses. We used the experimental data on ttA — ► J/ipX and 
fip — > ficcX to evaluate the parameters of the model and presented the relative contribu- 
tions of the photon-gluon fusion mechanism and the direct scattering on the constituent 
charmed quark (intrinsic charm) to the charmed structure function of the proton. 

We found that in irA collisions the so called " diffractive" component of J /if) can be 
well described by the coalescing of c and c constituent quarks. This success supports, 
in our view, Gribov's space-time picture of the hadron interaction with nuclei as well as 
the presence of long-living heavy quark fluctuations in the hadron ("intrinsic charm"). 
The results show also that the longitudinal distribution of the constituent heavy quarks 
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is harder than for light sea partons and has a shape like that of the valence quark distri- 
butions. 

From the comparison with the charm muon-pro duct ion at = 200 GeV we esti- 
mated the contribution from scattering on constituent charmed quarks to the total charm 
production. This contribution is about 1% and is expected to grow with the beam energy. 
At large values of the Bjorken variable x > 0.1 the scattering on the constituent charmed 
quark dominates the forward charm production in deep inelastic lepton-proton collisions. 

The considered approach taken for the excitation of Fock states, containing heavy 
quarks, can also be applied to states with beauty quarks in the hadron as well as to the 
hadronic component of the virtual photon ("resolved photon"). 

Finally we want to underline that the ep collider HERA is well suited to investigate 
heavy flavour production mechanisms over a wide kinematic region of (x, Q 2 ) which is 
inaccessible at other existing facilities. 
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APPENDIX 



The expression for the first-order radiative corrections <J^\z, A) to the charmed struc- 
ture function is: 

a ( i\z,\) = ^(5(l-2)|41nA-2 + v / TT4AL 

+ [3L 2 + 4L + 4L2 2 (^) + 2LlnA-4Lln(l + 4A) + 2L2 2 (f)]| 

+ 6Az(l -z)(3- 15z - 2z 2 + 8^ 3 ) + 4AV(8 - 77z + 65z 2 - 2z 3 ) + 
+ 16A 3 ^ 3 (1 - 21z + 12z 2 ) - 128AV] 



2L 



Vl+4Az 2 



[(1 + + 2^ 2 ) - 2A^(2 - llz - llz 2 ) - 8AV(1 - 9z)] 



8z 4 (l+A) 2 _ 4z 4 (l+2A)(l+4A) 2 L 
(1-2)+ Vl+4Az 2 (l- Z )+ 



(36) 



with the definitions: 



f = i i\z[l-{l-\)z\ , 

(1+2A2+V1+4A2 2 ) 2 ' 



v / l+4A+l . 
2 ' 



d = a-1; (37) 



In-- 



^2 = - /q 1 dz 
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Figure 1: Results of a fit to the NA3 data for J/ip production in up and pp interactions 
(a-e) using the model described in the text. The fit has been performed for nN collisions 
only, (f) The xq(x) distributions in the proton of the c (solid), valence u (dashed) quarks 
in the model and the PDF parameterization for valence u quark (dotted). 
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(b) 



Figure 2: Diagrams for photon-gluon fusion (a) and the scattering on the intrinsic 
charmed quark (b). 
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Figure 3: 0(a s ) corrections to the intrinsic charm structure function: (1), (2) 
bremsstrahlung and (3) — virtual gluon corrections. 
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Figure 4: Results of the fit of the sum PGF + 1C to the EMC data (a,b) and the ratio 
K = F 2 /C, /F 2 PGF of the contributions from intrinsic charm and photon gluon fusion (c). 
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